Itinerant Ferromagnetism in the Periodic Anderson Model 
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We introduce a novel mechanism for itinerant ferromagnetism, based on a simple two-band model. 
The model includes an uncorrelated and dispersive band hybridized with a second band which 
is narrow and correlated. The simplest Hamiltonian containing these ingredients is the Periodic 
Anderson Model (PAM). Using quantum Monte Carlo and analytical methods, we show that the 
PAM and an extension of it contain the new mechanism and exhibit a non-saturated ferromagnetic 
ground state in the intermediate valence regime. We propose that the mechanism, which does not 
assume an intra atomic Hund's coupling, is present in both the iron group and in some / electron 
compounds like Ce(Rhi_ 2 ,Ru a; )3B2 , LazCei-^RhaP^ and the uranium monochalcogenides US, USe, 
and UTe. 



I. INTRODUCTION 

Itinerant ferromagnetism was the first collective quan- 
tum phenomena considered as a manifestation of the 
strong Coulomb interactions which are present in an elec- 
tronic system. However, the origin of this phenomenon is 
still an open problem. Here we introduce a novel mech- 
anism for itinerant ferromagnetism, which is based on a 
simple two-band model. This work is an extension of 
a previous letter |jj where we have described the basic 
ideas. In this paper, we show that the mechanism is 
supported by the numerical results obtained from quan- 
tum Monte Carlo (QMC) simulations of the periodic An- 
derson model (PAM). We also analyse the experimen- 
tal consequences for some / electron compounds and the 
Iron group. However, before describing the details of our 
mechanism, it is useful to develop a historical perspective 
for the itinerant ferromagnetism. 

Seventy four years ago, Heisenberg 0] formulated his 
spin model to address the problem of ferroamgnetism, 
but as Bloch || pointed out, a model of localized spins 
cannot explain the metallic ferromagnetism observed in 
iron, cobalt and nickel. By introducing the effects of 
exchange terms into the Sommerfeld description of the 
free-electron gas, Bloch predicted a ferromagnetic ground 
state for a sufficiently low electron density. In 1934, 
Wigner 0] added the correlation terms to the model con- 
sidered by Bloch. After analyzing the effects of the cor- 
relations, he concluded that the conditions for ferromag- 
netism in the electron gas were so stringent as to be never 
satisfied. 

The first attempt at analyzing a real FM metal, like 
Ni, was made by Slater ||. He concluded that the main 
contribution to the exchange energy is provided by intra- 
atomic interactions. In the meantime, Stoner H intro- 



duced his picture where the metallic ferromagnetism re- 
sults from holes in the 3d band interacting via an ex- 
change energy proportional to the relative magnetization 
and obeying Fermi-Dirac statistics. However, the model 
considered first by Stoner [|| and later by Wohlfarth Q , 
did not take into account the correlations of the 3d elec- 
trons, except for the constraints imposed by the Pauli 
exclusion principle. In other words, they did not con- 
sider the fact that the Coulomb repulsion tends to keep 
the electrons apart. 

In 1953, the importance of these correlations was 
pointed out by van Vleck p5[ . He emphasized that the 
energy required to tear off an electron increases rapidly 
with the degree of ionization. (The energy of two Ni 
atoms in a 3d 9 configuration is appreciably lower than 
having one atom in the 3d 8 state and the other one in 
3d 10 .) Based on this observation, he proposed an alterna- 
tive picture (minimum polarity model) where the states 
of higher ionization in Ni are ruled out completely, and 
the configuration 3d 9 4 is considered to be 40 percent 3d 10 
and sixty percent 3d 9 . The lattice sites occupied by 3d 9 
and 3d 10 configurations are continuously redistributing 
in his picture. The van Vleck proposal is the precursor 
of the Hubbard model for infinite U. 

Following Slater (m), van Vleck |[5| speculated that 
the contamination by states of higher polarity, not in- 
cluded in his model, provides the exchange interaction 
(intra-atomic in this case) necessary for ferromagnetism. 
Hence he calculated the effective nearest-neighbor mag- 
netic interaction induced by second order perturbative 
fluctuations from the d 9 + d 9 configuration to the d 8 + d 10 . 
In this way, van Vleck arrived at a model which describes 
itinerant and correlated (only allowing d 9 and d 10 config- 
urations) holes with a nearest-neighbor exchange interac- 
tion aSi.Sj (generalized Heisenberg model). However, as 
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van Vleck explained at the end of his paper p5| , the sign 
and the magnitude of a are very sensitive to the precise 
values of the energies of the different possible intermedi- 
ate states (singlets or triplets) in the d 8 configuration. 

In 1963, the one-band Hubbard model was proposed 
independently by Gutzwiller |l6[ ], Hubbard JlTj and 
Kanamori |l8[ ] to explain the metallic ferromagnetism in 
the 3d transition metals. The Hubbard model incorpo- 
rates the kinetic energy in a single nondegenerate band 
with an intra-atomic Coulomb repulsion U to describe 
the electrons in the s band of the transition metals. In 
contrast to the previous models, the Hubbard model does 
not include any explicit exchange interaction which favors 
a ferromagnetic phase. The implicit question raised by 
this proposal is: Can ferromagnetism emerge from the 
interplay between the kinetic energy and the Coulomb 
repulsion, or it is strictly neccesary to include an ex- 
plicit exchange interaction provided by the intra-atomic 
Hund's coupling? This simple question becomes even 
more relevant if we consider /-electron itinerant ferro- 
magnets, like CeRli3B2 fl9fl , whose only local magnetic 
coupling is antiferromagnetic. 

Unfortunately, with the exception of Nagaoka's |2(]] 
and Lieb's pljp2j theorems, the subsequent theoreti- 
cal approches were not controlled enough to determine 
whether the Hubbard model has a FM phase. The cen- 
tral issue is the precise evaluation of the energy for the 
paramagnetic (PM) phase. Because it does not properly 
incorporate the correlations, mean field theory overes- 
timates this energy and predicts a large FM region p3| . 
In contrast, numerical calculations have narrowed the ex- 
tent of this phase to a small region around the Nagaoka 
point @. 

There are several controversies in the history of itiner- 
ant ferromagnetism. Most of them originate in the lack of 
reliable methodologies to solve strongly correlated Hamil- 
tonians. In addition, the number of compounds exhibit- 
ing ferromagnetism has increased during the last several 
decades. There is a new list of Ce and U based com- 
pounds which are itinerant ferromagnets [^) with large 
Curie temperatures (of the order of 100°K) and the only 
explicit magnetic interaction is antiferromagnetic (Kondo 
like). This observation stimulated us to reconsider the 
situation of the iron group and ask whether the ferromag- 
netism originates in the intra-atomic Hund's interaction 
or is just a consequence of the strong Coulomb repul- 
sion between electrons with a particular band structure. 
One possible and plausible alternative is that these two 
phenomena are cooperating to stabilize a ferromagnetic 
ground state. In any case, it is important to elucidate 
whether a realistic lattice model only containing repul- 
sive terms (for instance, the Hubbard model) can sustain 
a ferromagnetic state, or it is necessary to invoke ad- 
ditional terms which are explicitly ferromagnetic. The 
Hubbard Hamiltonian is in the first category. However, 
the most accurate numerical calculations seem to indicate 



that this model does not sustain a ferromagnetic ground 
state. 

Going beyond the simple one-band Hubbard model 
has been advocated, for instance, by Vollhardt et al 
p3[ . They note that the inclusion of additional Coulomb 
density-density interactions, correlated hoppings, and di- 
rect exchange interactions favors FM ordering in the 
single-band Hubbard model. In fact, a very simple analy- 
sis shows that increasing the density of states T>(E) below 
the Fermi energy Ep and placing Ep close to the lower 
band edge increases the FM tendency. One can achieve 
this by including a next nearest neighbor hopping t' or 
by placing the hoppings on frustrated (non-bipartite) lat- 
tices. The effectiveness of t' was studied numerically by 
Hlubina et. al. J2^] for the Hubbard model on a square 
lattice. They found a FM state when the van Hove sin- 
gularity in T>(E) occurred at Ep. However, this phase 
was not robust against very small changes in t' . 

After seven decades of intense effort, the microscopic 
mechanisms driving the metallic FM phase are still un- 
known |^3|,^8| . We still do not know what is the minimal 
lattice model of itinerant ferromagnetism and, more im- 
portantly, the basic mechanism of ordering. 

While the Hubbard model is so reluctant to have a FM 
state, there is an increasing amount of evidence indicat- 
ing that the Periodic Anderson Model (PAM) has a FM 
phase in a large region of its quantum phase diagram 
p],E9-39 . Since the d orbitals of the transition metals 



are hybridized with the s — p bands, we can consider the 
inclusion of a second band as the next step in the search 
of itinerant ferromagnetism from pure Coulomb repul- 
sions. Indeed a very simple extension of the PAM can 
be used to describe the physics of the iron group |^0|,[30| . 
In addition, there is large number of cerium and ura- 
nium compounds, like Ce(Rhi_ x Ru I )3B2 (T c = 115K) 
HH, CeRh 3 B 2 ||, and US, USe, and UTe @ which 
are metallic ferromagnets and can be describred with the 
PAM. 

Ferromagnetism is readily found in the PAM by mean- 
field approximations in any dimensions p3|-^9|. Using a 
slave-boson mean-field theory (SBMFT) for the symmet- 
ric PAM, Moller and Wolfe (33) found a PM or antifer- 
romagnetic (AF) phase at half filling depending on the 
value of the Coulomb repulsion U. By lowering the den- 
sity of electrons from 1/2 filling, they also found a smooth 
crossover from AF to FM order via a spiral phase. Just 
before 1/4 filling, they got a first-order transition from 
FM to AF order. More recently, the SBMFT calculations 
of Doradzihski and Spalek |54],[35| found wide regions of 
ferromagnetism in the intermediate valence regime that 
surprisingly extended well below 1/4 filling. 

A ferromagnetic phase is also obtained when the dy- 
namical mean-field theory (DMFT) is applied to the 
PAM §§-|3j|. Tahvildar-Zadeh et al. found a region of 
ferromagnetism and studied its temperature dependence. 
At very low temperatures, their ferromagnetic region ex- 
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tended over a wide range of electron fillings and in many 
cases embraced the electron filling of 3/8. They pro- 
posed a specific Kondo-induced mechanism for ferromag- 
netism at 3/8 filling that has the conduction electrons 
in a spin-polarized charge density-wave anti-aligned with 
the ferromagnetically aligned local moments on the va- 
lence orbitals. More recently Meyer and Nolting 37 3H| ] 
appended perturbation theory to DMFT and also pre- 
dicted ferromagnetism over a broad range of electron fill- 
ing extending below 1/4 filling. In addition Schwieger 
and Nolting also considered an extension of the PAM, 
similar to the one considerd here, to estimate the impor- 
tance of s — d hybridization for the magnetic properties 
of transition metals. 

There is also a considerable amount of numerical ev- 
idence showing ferromagnetic solutions for the ground 
state of the PAM. Noack and Guerrero pi]] , for example, 
found partially and completely saturated ferromagnetism 
using the density matrix renormalization group (DMRG) 
method in one dimension. They considered a parameter 
regime where there is one electron in each / orbital. For a 
sufficiently large value of U, the model exhibited a ferro- 
magnetic ground state. Beyond an interaction-dependent 
value of the doping and a doping-dependent value of U, 
this state disappeared. The ferromagnetic phase was a 
peninsula in a phase diagram that was otherwise a sea of 
paramagnetism except at 1/4 and 1/2 filling where the 
ground state of the PAM was antiferromagnetic. 

Our previous [^9| and new QMC results qualitatively 
agree with the DMRG work; however, the phases we 
find quantitatively and qualitatively disagree with those 
derived from the mean-field approximations. Quantita- 
tively, we find ferromagnetism in a narrower doping range 
than the one predicted by the DMFT and SBMFT calcu- 
lations. For fillings between 3/8 and 1/2, QMC predicts 
a PM region, whereas mean-field theory predicts ferro- 
magnetic states in part of that region. In fact, at a filling 
of 3/8 where DMFT calculations predict ferromagnetism, 
we find a novel ground state of an entirely different sym- 
metry. Instead of ferromagnetism, QMC finds a resonat- 
ing spin density-wave (RSDW) state; that is, the ground 
state was a linear combination of two degenerate spin- 
density waves characterized by the (-7T, 0) and (0, tt) wave 
vectors. 

The novel mechanism we introduce in the present pa- 
per operates when the system is in a mixed valence 
regime. This regime has been studied numerically only 
in the context of DMFT |j7|. We will show however that 
the ferromagnetic solution obtained with DMFT in the 
mixed valence regime has a different origin and therefore 
is not representative of our new mechanism. The main 
ingredients for our mechanism are an uncorrelated dis- 
persive band which is hybridized with a correlated and 
narrow band. We show the PAM supports our mecha- 
nism by doing quantum Monte Carlo (QMC) simulations 
on one and two dimensional lattices. The results of these 



simulations are interpreted with an effective Hamiltonian 
derived from the PAM. In this way, we establish that the 
new mechanism can be interpreted as a generalization to 
the lattice of the first Hund's rule for the atom. The two 
level band structure generated by the gap recreates for 
the lattice, the shell-like level strucure of the hydrogenic 
atom. When the lower shell is incomplete, the local part 
of the Coulomb interaction is minimized by polarizing 
the electrons which are occupying the incomplete shell. 



II. MODEL 

The PAM was originally introduced to explain the 
properties of the the rare-earth and actinide metallic 
compounds including the so called heavy fermion com- 
pounds. A very simple extension of this model can also 
be applied to the description of many transition metals 
p6| , ^0| . The basic ingredients of this model are a narrow 
and correlated / band hybridized with a despersive and 
uncorrelated d band. The Hamiltonian associated with 
this model is: 
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where b\ a and a\ a create an electron with spin a in b 
and a orbitals at lattice site r and n" CT = a\ a a Ta - The tb 
and t a hoppings are only to nearest-neighbor sites. When 
t a = 0, the Hamiltonian is the standard PAM. For the 
/ electron compounds, the a and b orbitals play the role 
of the / and d orbitals, and t a « 0. For transition met- 
als, they correspond to the 3d and 4s orbitals. Unless 
otherwise specified, we will set tb = 1. 

For [7 = 0, the resulting Hamiltonian Hq is easily di- 
agonalizcd: 



(2) 



where the dispersion relations for the upper and the lower 
bands are: 



1 



+ e a k ±y/(el-eZ)*+4Vi 



(3) 



3 



with 
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e£ = -2t b y^cos k Xi , 

i=l 

D 

e k = e a ~ 2t Q ^2 C0S k * 
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for a hypercubic lattice in dimension D. The operators 
which create quasi-particles in the lower and upper bands 
are: 



(5) 



with 
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The noninteracting bands £^ are plotted in Fig. 1 for a 
one dimensional system. If |V| -C |t&|, we can identify 
regions with well defined / or d character in the lower 
and the upper bands. In particular, the case illustrated 
in Fig. 1 corresponds to a situation where the a and the b 
bands were crossing before being hybridized (e a > — 2tb). 
We see a small region in the center of the lower band 
which is dispersive and large regions on both sides which 
are nearly flat. The upper band exhibits the opposite 
behavior. The nearly flat regions in each band corre- 
spond to states with a predominant / character, while 
the dispersive regions are associated to the states with d 
character. 




FIG. 1. Illustration of the effective model and the FM 
mechanism. A is the hibridization gap and 8 a is the inter- 
val of energy where the electrons are polarized. 



III. NEW MECHANISM FOR 
FERROMAGNETISM 

The PAM has different regimes depending on the val- 
ues of its parameters and the particle concentration 
n = N e /AN (N e is the total number of particles and N is 



the number of unit cells). If V <C |e a | and V <C \U + e a \, 
there is one particle (magnetic moment) localized in each 
a orbital, and the fluctuations to the conduction band 
can be considered in a perturbative way. By this pro- 
cedure, the PAM can be reduced to the Kondo Lattice 
Model (KLM) H] which contains only one parameter 
Jx/t (with Jk = and t = if,). The KLM has been exten- 
sively studied ^ , and the evolution of its phase 
diagram is described for instance in a review article by 
Tsunetsugu et al [Q . One of the earliest approaches to 
the KLM is the mean field treatment of Doniach pj| for 
the related one dimensional Kondo necklace. For half fill- 
ing, this approximation leads to a transition from a Neel 
ordered state in the weak coupling regime (Jk <C \t\) 
to a nonmagnetic 'Kondo singlet' state above the critical 
value J K = t. 

Lacroix and Cyrot [fl9| did a more extensive mean field 
treatment for three dimensional KLM. They also found 
a magnetically ordered state for weak coupling. For low 
density of conduction electrons. In their phase diagram, 
the ordered state is ferromagnetic for low and intermedi- 
ate densities of conduction electrons, and antiferromag- 
netic in the vicinity of half filling. The 'Kondo singlet' 
phase appears above some critical value of Jj^i 71 ) in the 
whole range of concentrations. 

Using another mean field treatment for the one dimen- 
sional KLM, Fazekas and Miiller-Hartmann |5(J obtained 
a phase diagram containing only magnetically ordered 
phases: spiral below some critical value of Jx/t which 
depends on the particle density and ferromagnetic above 
this value. To get this result, they fixed the orientation 
of the localized spins in a spiral ordering and minimized 
the total energy with respect to the wave vector of the 
spiral. Even though this treatment of the spin polarized 
state is valid for classical spins, it neglects completely 
the Kondo singlet formation which occurs in the strong 
coupling limit for the considered case (S = 1/2). 

Sigrist et al |5l|| gave an exact treatment of one di- 
mensional KLM for the strong coupling regime Jk 3> t 
finding a ferromagnetic phase for any particle density. 
However, it is important to remark that the mechanism 
driving the ferromagnetism in the later case is not the 
same as the double-exchange mechanism associated with 
the mean field solution of Fazekas and Miiller-Hartmann 
|p0| . To understand this difference, we just need to notice 
that for Jif/i = oo mean field Mj predicts a ferromag- 
netic solution while the exact solution has a complete 
spin degeneracy. Therefore double-exchange is not the 
mechanism driving the ferromagnetic phase of the KLM 
(at least when the localized spins are S = 1/2). 

The real mechanism has been unveiled by Sigrist et al 
[ mj who used degenerate perturbation theory to deter- 
mine the lifting of this degeneracy when the ratio Jif/t 
becomes finite. The new ground state is an itinerant fer- 
romagnet for any concentration of conduction electrons. 
In this state, the spins which are not participating in 
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Kondo singlet are fully polarized. We can see from their 
solution that the motion of the Kondo singlets stabilizes 
the FM state in a way similar to Nagaoka's solution pOl . 
The second order effective Hamiltonian obtained after the 
perturbative calculation includes nearest-neighbor hop- 
ping t/2, plus a next-nearest-neighbor correlated hopping 
t' which is order t 2 / Jk- Then there are two different ways 
to move a Kondo singlet from one site to its next-nearest- 
neighbor: by two applications of t/2 or by one applica- 
tion of t' . Only when the background is ferromagnetic do 
both processes lead to the same final state. If t' has the 
appropriate sign (which is the case for the KLM |5l| ) the 
resulting interference is "constructive" and the FM state 
has the lowest energy. We can see in this example that 
the motion of the Kondo singlet can stabilize a magnetic 
phase. 

There arc different regimes for which the PAM cannot 
be reduced to a KLM by a perturbative approach. One of 
these situations corresponds to the intermediate valence 
region: e a ~ Ep. In this case nf is not longer close to 
one and the a electrons can move. 

In Fig. I, we illustrate the (one-dimensional) non- 
interacting bands for the case of interest: e a close to Ep 
and above the bottom of the b band. If | | <§; |tf,|, we 
can identify two subspaces in each band where the states 
have either b (ip subspace) or a (<p subspace) character. 
The size of the crossover region around the points where 
the original unhybridized b and a bands crossed is pro- 
portional to |Vyt&|; that is, it is very small. The creation 
operators for the Wannier orbitals tp T(T and (j> T(T associated 
with each subspace are: 
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where N is the number of sites. The subsets K> and 
K< are defined by: K> = {k : |uk| > |«k|} and 
K< = {k : |vk| > |uk}- Since the ip and the (j> subspaces 
are generated by eigenstates of H , it is clear that both 
subspaces can only be mixed by the interacting term Hjj. 
Therefore in the new basis we have: 

Ho = H+ + H* = E rtAU^a + E itAl^ 



The segmented structure of the ip and the <f> bands intro- 
duce oscillations in the hoppings rf and rf as a function 
of the distance |r|. 

Because the U term in H involves only the a orbitals, 
the matrix elements of H connecting the <f> and ip sub- 
spaces are small compared to the characteristic energy 
scales of the problem (the matrix elements of H within 
the subspaces). To see this we express a\ a as a function 
of (f)\. a and V4<7 by first inverting Eqs. |^ and |?j to find: 



(10) 



where the weights W r and w T are defined by: 
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The value of these weights as a function of the distance 
r is plotted in Fig. 2 (for e a = —1.5, V = O.lt and 
t a = 0). From Fig. 2 we can see that Wo is much larger 
than any other weight. This is so because the cj> orbitals 
have predominantly a character, while the tp orbitals have 
mostly b character in the considered region of parameters. 
Therefore we can approximate the creation operator a\ a 
by: 



E 



(12) 



As a consequence of this approximation, the a subspace 
becomes invariant under the application of H. In ad- 
dition, because \W \ 3> |W r j | (see Fig. 2), we can es- 
tablish a hierarchy of terms where the lowest order one 
corresponds to a simple on-site repulsion: 



with the hoppings and rf given by the following ex- 
pressions: 
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with U = C/|IF | 4 and 



r,(7rr.fT " 



The next order 



terms, containing three and two Wo factors, are much 
smaller and are essentially the same as the intersite in- 
teractions which in the past were added to the Hubbard 
model to enhance the ferromagnetism p3[. 
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FIG. 2. Weights W T and w r as a function of the distance 
|r| for a one dimensional system. 

Adding to Hq we get the effective Hamiltonian: 

r,r' ,a 

+ uY,4Ai ( 14 ) 

r 

The ip and cj> orbitals form uncorrelated and correlated 
non-hybridized bands: H e ff = + H^. For the 4> or- 
bitals we obtain an effective one band Hubbard model 
with the peculiar double shell like dispersion relation 
shown by the thick lines in Fig. 1. 

Particularly for t a — 0, has a very large density of 
states in the lower shell of the <p band J2^] which is located 
near e a . From Fig. 1 it is also clear that the electrons 
first doubly occupy the uncorrelated ip band states which 
are below e a . However, when Ep gets close to e a , i.e., the 
system is in the mixed valence regime, the electrons close 
to the Fermi level go into some of the correlated <fi states. 
Then, the interaction term j , combined with the dou- 
ble shell band structure of Hq , gives rise to a FM ground 
state (GS): The electrons close to Ep spread to higher 
unoccupied k states and polarize, which causes the spa- 
tial part of their wave function to become antisymmetric, 
eliminating double occupancy in real space and reducing 
the Coulomb repulsion to zero. The cost of polarizing 
is just an increase in the kinetic energy proportional to 
5 a r-j hvpSk, where vp is the Fermi velocity and 5k is the 
interval in k space in which the electrons are polarized. 

To determine the stability of this unsaturated FM 
state, we compare its energy with that of the PM state. 
If we were to build a nonmagnetic state with only the 
states of the lower (f> shell, we would find a restricted de- 
localization for each electron because of the exclusion of 
the finite set of band states (k-states) in the upper shell. 
To avoid the Coulomb repulsion U for double occupying 
a given site, the electrons need to occupy all k-states. 



This means they have to occupy the 4> states in the up- 
per and lower shells. This restricted derealization is a 
direct consequence of Heisenberg's uncertainty principle, 
and the resulting localization length depends on the wave 
vectors, where the original b and a bands crossed, that 
define the size (Ak) of each shell. The energy cost for 
occupying the 4> states in the upper shell is proportional 
to the hybridization gap A. Therefore if U is the dom- 
inant energy scale in the problem and A ^S> S a , the FM 
state lies lower in energy than the nonmagnetic state. 
Under these conditions, the effective FM interaction is 
proportional to the hybridization gap A. 

This mechanism for ferromagnetism on a lattice is 
analogous to the intra-atomic Hund's mechanism polar- 
izing electrons in atoms. In atoms, we also have differ- 
ent degenerate (the equivalent of S a is zero) shells sepa- 
rated by an energy gap. If the valence shell is open, the 
electrons polarize to avoid the short range part of the 
Coulomb repulsion (again reflecting the Pauli exclusion 
principle) . The energy of an eventual nonmagnetic state 
is proportional either to the magnitude of the Coulomb 
repulsion or to the energy gap between different shells. 
The interplay between both energies sets the scale of 
Hund's intra-atomic exchange coupling. 

IV. NUMERICAL METHOD 

Our numerical method, the constraincd-path Monte 
Carlo (CPMC) method, is extensively described and 
benchmarked elsewhere. |5^^] Here we only discuss its 
basic features, assumptions and special details about our 
use of it. 

In the CPMC method, the ground-state wave function 
\^o) is projected from a known initial wave function \^t) 
by a branching random walk in an over-complete space 
of Slater determinants In such a space, we can write 
l^o) = S^c^l^)' where > 0. The random walk pro- 
duces an ensemble of \4>), called random walkers, which 
represent \^o) in the sense that their distribution is a 
Monte Carlo sampling of c^/^^c^, that is, a sampling 
of the ground-state wave function. 

To completely specify |\&o)j only determinants satisfy- 
ing (^>o\(f>) > are needed because \^q) resides in either 
of two degenerate halves of the Slater determinant space, 
separated by a nodal plane. In the CPMC method the 
fermion sign problem occurs because walkers can cross 
this plane as their orbitals evolve continuously in the ran- 
dom walk. Without a priori knowledge of this plane, we 
use a trial wave function \^t) and require (^p\<l>) > 0. 
The random walk solves Schrodinger's equation in de- 
terminant space, but under an approximate boundary- 
condition. This is what is called the constrained-path 
approximation. 

The quality of the calculation depends on the quality 
of the trial wave function \^t)- Fortunately, extensive 
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testing has demonstrated a significant insensitivity of the 
results to reasonable choices: Since the constraint only 
involves the overall sign of its overlap with any deter- 
minant some insensitivity of the results to \^t) is 
expected g|^6lj. 

Besides as starting point and as a condition constrain- 
ing a random walker, we also use \^t) as an importance 
function. Specifically, we use (^>t\4>) to bias the ran- 
dom walk into those parts of Slater determinant space 
that have a large overlap with the trial state. For all 
three uses of \^t), it clearly is advantageous to have 
I^t) approximate \^q) as closely as possible. Only in 
the constraining of the path does \^t) 7^ l^o) in general 
generate an approximation. 

We constructed |?/>t) = Y\ a \4>t) from the eigenstates 
of the non-interacting problem. Because the total and z- 
component spin angular momentum, S and S z , are good 
quantum numbers, we could choose unequal numbers of 
up and down electrons to produce trial states and hence 
ground states with S = S z = |(iV| — N{). Whenever 
possible, we would simulate closed shells of up and down 
electrons, as such cases usually provided energy estimates 
with the least statistical error, but because we wanted to 
study the ground state energy as a function of S, we 
frequently had to settle for just the up or down shell be- 
ing closed. In some cases, the desired value of S could 
not be generated from either shell being closed. Also 
we would select the non- interacting states so \ipr) would 
be translationally invariant, even if these states used did 
not all come from the Fermi sea. The use of unrestricted 
Hartree-Fock eigenstates to generate|0y) instead of the 
non-interacting eigenstates generally produced no signif- 
icant improvement in the results. 

In particular, we represented the trial wavefunction as 
a single Slater determinant whose columns are the A CT 
single-particle orbitals obtained from the exact solution 
of Hq. We chose the orbitals with lowest energies given 
by and filled them up to a desired number of electrons 
N e . 

i^>=rKjo>> as) 

k,cr 

where |0) represents a vacuum for electrons. Since our 
calculations were performed for a less than full lower 
band, only states from the lower band were used to con- 
struct the trial wavefunction. 

In a typical run we set the average number of ran- 
dom walkers to 400. We performed 2000 Monte Carlo 
sweeps before we taking measurements, and we made the 
measurements in 40 blocks of 400 steps. By choosing 
At = 0.05, we reduced the systematic error associated 
with the Trotter approximation to be smaller than the 
statistical error. In measuring correlation functions, we 
performed between 20 to 40 back-propagation steps. 



V. QUANTUM MONTE CARLO RESULTS 

In section III, we have described a new mechanism for 
itinerant ferromagnetism which is present in the mixed 
valence regime for n > 1/4. In addition, we mentioned 
that the system is also expected to be FM when the / 
magnetic moments are localized (\V\ <C \Ep — e a \) be- 
cause the effective RKKY coupling is negative when the 
Fermi surface is small {kp ~ 0). In this section we show 
that the itinerant and the localized ferromagnetic states 
are continuously connected in the phase diagram of the 
PAM. However, the energy scale of the first state is much 
larger than the RKKY interaction which characterizes 
the second one. The existence of a crossover region be- 
tween both states could explain the fact that there are 
some /-electron compounds for which it is very difficult 
to determine whether they are itinerant or localized fer- 
romagnets. 

In addition, we will see that the QMC results are con- 
sistent with the simple picture derived from our effective 
model (Eq. Hj). According to that picture, the ferro- 
magnetic state in the mixed valence regime should be 
similar to a partially polarized non-interacting solution 
where the polarized electrons are the ones occupying the 
/-character orbitals. It is only in the crossover region of 
size V 2 /tb (see Fig. |l|), where the orbitals have a mixed 
character, that the correlations introduce an appreciable 
effect. This effect is the well known Kondo-like singlet 
correlation between the d and the / electrons. How- 
ever, it is important to remark that these Kondo sin- 
glets only exist in an energy interval V 2 /W (where W 
is the d bandwidth), and therefore the number of Kondo 
singlets Nks is much smaller than the number of mag- 
netic moments N M m- N ks /N M m ~ V 2 /W 2 . This is 
a simple manifestation of the "exhaustion" phenomenon 
described by Nozieres |32|,^3| . Since most of the / mag- 
netic moments are ferromagnetically polarized, the role of 
these few Kondo singlets is marginal in our FM solution. 
Therefore, for the mixed valence regime with n > 1/4, 
the / magnetic moments which are not screened by d 
electrons develop an effective magnetic interaction as a 
consequence of the interplay between the local Coulomb 
interaction and the particular band structure. In other 
words, the "collective Kondo state" which was proposed 
in the past ]62| , |63[ | is replaced by band ferromagnetism 

In fact, the nature of local moment compensation in 
the PAM differs qualitatively from that in the single im- 
purity Anderson model (^9|. In the PAM, if the ground 
state is a singlet, then 

( s f(tf) = -E< 5 scw)) - EWW)) 

(16) 

In the impurity model, the last term is absent, and the 
resulting expression is the analytic statement of the well- 



7 



known Clogston-Anderson compensation theorem that 
express the compensation of the f-moment by the con- 
duction electrons. In the PAM, on the other hand, the 
last term dominates the first so the f-moment is compen- 
sated largely by correlations induced among themselves. 

To understand the nature of the FM solution, we plot- 
ted the mean occupation number of the quasiparticle 
operators which diagonalizc the non-interacting problem 
(see Eq. ||) . This is shown in Fig. || for the lowest energy 
PM (S = 0) state in a two dimensional cluster of 12 x 12 
unit cells. There we can see that the d-character states of 
the lower band are close to be doubly occupied. In con- 
trast, the /-character region has a smaller occupation 
number (~ 0.7). It is remarkable that the populations of 
the /-character states in the lower and the upper bands 
are very similar. This derealization in the momentum 
space is a direct consequence of the tendency to avoid 
double-occupancy in the real space. This tendency indi- 
cates that the energy increase of the PM state due to the 
inclusion of U is proportional to the hybridization gap. 



1.5 




FIG. 3. Mean value of the occupation numbers of the 
non-interacting band states for the paramagnetic S — solu- 
tion of the PAM in a two dimensional lattice. Legends: n;„ (k) 
represents lower-band occupation number of up-spins, n„<j(k) 
represents upper-band occupation number of down-spins, etc. 

Fig. ^ shows the quasiparticle occupation numbers for 
the partially saturated FM ground state. While the b- 
character states of the lower band are still close to be- 
ing doubly occupied, the a-character ones are polarized 
and the occupation number is one for any of them. In 
contrast to the PM solution (see Fig.[|), the occupation 
number for the /-character states in the upper band is 
much smaller. This difference can be understood in the 
following way: the polarized /-electrons can localize in 
momentum space because the Pauli exclusion principle 
prevents double occupancy in real space (the spatial wave 
function is completely antisymmetric). Therefore the a- 
electrons do not need to occupy the upper band states 
and the energy increase due to the repulsive U term is 
proportional to hvF&k- The non-zero amount of electrons 



occupying the center of the upper band (see Fig. [|) is re- 
lated to the crossover region for which the a and the b 
character of the states are comparable. Since the elec- 
trons occupying these states are not polarized, the effect 
of the Coulomb repulsion U is the transfer of spectral 
weight from the lower to the upper band to avoid double 
occupancy in real space. 

We can see from Figs. || and ^ that a similar behavior is 
obtained for one dimensional systems. As we explained in 
section III, the mechanism for this ferromagnetism works 
in any finite dimension. Notice in Fig. |^ that there is a 
jump in the occupation number of the lower band as a 
function of k. The inverse of this jump is proportional 
to effective mass of the quasiparticles of the pramagnctic 
solution. When e a decreases, the system evolves into a 
state where the a electrons are localized. This evolution 
is reflected in the decrease of the jump and the conse- 
quent increase of the effective mass. 
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FIG. 4. Mean value of the occupation numbers of the 
non-interacting band states for the ground state (partially 
saturated ferromagnet with M = S/N = 0.43) of the PAM in 
a two dimensional lattice. 
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FIG. 5. Mean value of the occupation numbers of the 
non-interacting band states for the paramagnetic S — solu- 
tion of the PAM in a chain. 
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FIG. 6. Mean value of the occupation numbers of the 
non-interacting band states for the ground state (partially 
saturated ferromagnet with M = S/N = 0.3) of the PAM in 
a chain. 



For U = the system is a PM metal. Therefore there 
must be a critical value of the on site repulsion U c which 
separates the PM from the FM region. The value ob- 
tained for U c is: U c ~ 0.25 ~ A = 0.242 for V = 0.5 
and U c ~ 0.11 ~ A = 0.0981 for V = 0.31 (see Fig. 0). 
This is also in agreement with the mechanism described 
in section III. If A ^> 6 and U becomes larger than the 
hybridization gap A, the system evolves into a FM state 
to avoid the double occupancy without an increase in the 
kinetic energy proportional to the hybridization gap. 



FIG. 7. a) Energy difference between partially polarized 
FM ground state and the lowest energy paramagnetic state 
as a function of U. b) Magnetization as a function of U. 
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FIG. 8. Minimum energy as a function of the magnetiza- 
tion. 



According to Figs, ^b and || the magnetization seems to 
increase gradually when U is increased beyond its critical 
value. This behavior suggests that the FM transition as 
function of U is of second order. If this is so, the proper- 
ties of the PM Fermi liquid which is obtained for U < U c 
should be strongly affected by the FM fluctuations when 
U approaches U c . It is known that the effective mass 
of the quasi-particles diverges in the approach to a zero 
temperature ferromagnetic instability J52| . In other 
words, the behavior of the PM Fermi liquid cannot be un- 
derstood by analogy with the one impurity problem. The 
Kondo temperature, which is the characteristic energy 
scale of the one impurity problem, is replaced by a new 
Fermi temperature which is dominated by the ferromag- 
netic fluctuations and goes to zero when U approaches 
to U c from below. 



Another relevant parameter for the FM solution is the 
hybridization V. For V = and e Q = Ep there is a 
complete spin degeneracy for the electrons occupying the 
localized a orbitals. By increasing V, we are simultane- 
ously changing the Fermi velocity vp and the hybridiza- 
tion gap A. For small values of V, A is much larger 
than S — hvpSk in the region under consideration. For 
this reason, a non-zero value of V removes the original 
spin degeneracy stabilizing the partially polarized FM 
solution (see Fig. |J). When V is larger than t, the two 
relevant energy scales, A and 5, become of the same or- 
der and the partially polarized FM is replaced by a PM 
phase. In the unrealistic large V limit (\V\ ^> \tf>\, U, |e a |), 
the ground state consists of local Kondo singlets moving 
in a background of localized spins (5(| . 
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FIG. 9. a) Energy difference between the FM ground state 
and the lowest energy paramagnetic state as a function of V. 
b) Magnetization as a function of V. 

Finally, the most sensitive parameter for the stabi- 
lization of the FM state is the difference |e„ — Ep\. In 
Fig. [ll], we show the energy difference \AE\ between the 
FM ground state and the lowest energy PM (S = 0) state 
as a function of e a . When e a is considerably smaller than 
Ef, the a electrons are localized and the magnetism is 
dominated by the RKKY interaction. The order of this 
interaction is V 4 . This gives the small value of \AE\ 
when the a levels are below the bottom of the conduc- 
tion band e a < — 2D\tb\. The most stable region for 
the FM state (maximum value of |Ai?|) starts when e a 
reaches the Fermi level. For the case of Fig. [lC], this oc- 
curs at e a ~ — 1.9t. Again this result is in agreement 
with the mechanism described in section III. If we con- 
tinue increasing the value of e a , the number of a electrons 
decreases and the magnetization is consequently reduced 
(see Fig. [Tojo) . Finally, when e is no longer close to the 
bottom of the conduction band, 8 a becomes comparable 
to A and the ferromagnetism disappears. 

Fig. |ll] shows the e a dependence of \AE\ for two and 
three dimensional clusters. As in the one dimensional 
case, the stablity of the FM phase increases when the sys- 
tem appraches the mixed valence regime. If e a is further 
increased, the magnetization goes to zero in a smooth 
way indicating that the associated quantum phase tran- 
sition is of second order. 

In Figs. [[2] and [l3| we show the scaling of AE and the 
magnetization per site in one and two dimensional sys- 
tems. The 1/N extrapolation for the one dimensional 
systems indicates that the FM ground state is stable in 
the thermodynamic limit. The extrapolated value for 
the magnetization is ~ 0.3, i.e., only a fraction of the 
/ electrons is polarized. The size effects are stronger in 
two dimensional systems. In Fig. we show two differ- 
ent cases: the circles correspond to a sequence of tilted 
suqare clusters which allows us to fix the concentration in 
n = 0.3; the squares correspond to a sequence of untilted 



square clusters for which the value of the concentration is 
the closest to n = 0.29. Despite the considerable size ef- 
fects, these results indicate that the FM state is stable in 
the thermodynamic limit. In this case, the extrapolated 
magnetization is close to 0.4. 
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FIG. 10. a) Energy difference between FM ground state 
and the lowest energy paramagnetic state as a function of e a 
for a one dimensional system, b) Magnetization as a function 
of e a . 
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FIG. 11. a) Energy difference between FM ground state 
and the lowest energy paramagnetic state as a function of e a 
for two and three dimensional systems, b) Magnetization as 
a function of e a - 

Meyer and Nolting have also found a FM solu- 
tion for the PAM in a similar region of parameters us- 
ing DMFT. However it is important to remark that the 
mechanism for ferromagnetism described in section III 
only works in finite dimension because the volume of the 
upper shell relative to the volume of the lower one goes to 
zero when the dimension goes to infinity. Therefore, the 
reason why a FM state is obtained with DMFT should 
be different. This is also reflected by the fact that the en- 
ergy scale (T c ) of the FM solution found with DMFT for 
the localized regime is larger than the one for the mixed 
valence regime (region IV of Ref. j37J). This behavior is 
opposite to our result (see Fig. [To| ). 
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VI. EXPERIMENTAL CONSEQUENCES 
A. Cerium Compounds 

During the last few years new experimental results 
have confirmed that there are Ce based compounds which 
cannot be treated as typical Kondo systems. For in- 
stance, CeRh3B2 has a very high FM ordering temper- 
ature (T c = \\h°K) which sicks out from a localized 
4/-electron description [42| . In addition the alloy series 
CetRh^RoJaBa @ and La^Cei.^RhgBa @ exhibit 
many unusual characteristics which require a new macro- 
scopic description with respect to the competition among 
classical Kondo versus RKKY interactions |^5| . 

Absorption edge spectroscopy measurements of 
Ce(Rhi_ K Ru x )3B2 for different values of x indicate that 



the stoichiometric compound CeRli3B2 is in the mixed 
valence regime (fluctuating between the 4/ 1 and the 
4/° configurations) |Q. After doping with Ru, there is 
strong transfer of weight from the 4/ 1 line to the 4/° 
structure. This change can be understood in the context 
of the PAM if we take into account that the volume of 
the system decreases when Ce is replaced by Rh Q . In 
this situation the width of the conduction band increases 
and some / electrons are transferred to d character or- 
bitals (see Fig.^). According to our results this change 
must decrease the value of the zero temperature mag- 
netization and the Curie temperature T c . By increasing 
the doping level we can reach a situation where most of 
the / electrons that were polarized in the stoichiometric 
compound are transferred to the d character orbitals and 
the zero temperature magnetization is very small. In 
this limit the system should have some PM states very 
close to the Fermi level (see Fig. |l|). When the energy 
difference between the Fermi level and these PM states 
becomes smaller than T c , the magnetization can increase 
with temperature because the electrons which are occu- 
pying the PM states near the Fermi level are thermally 
promoted to the / character states which are above the 
Fermi level. In this process, the electrons are polarized 
because of the mechanism discussed above. This ex- 
plains the finite temperature peak in the magnetization 
of CetRhi^Ru^Ba (for x between 0.06 and 0.125) 
pj"| , [i"2"| that suggests an ordered state with high entropy. 
The source of the large entropy is thus associated with 
charge and not with spin degrees of freedom which is 
why a state with larger M has a higher entropy. From 
this analysis we predict that the integral of the entropy 
below T c , which can be extracted from the specific heat 
measurements, contains a considerable contribution from 
the the charge degrees of freedom. 

When Ce is replaced by La, the volume of the system 
increases |6~i|l and the magnetic moments become more 
localized. In this case, the weight in the absorption edge 
spectroscopy is transferred from the 4/° structure to the 
4/ 1 line. Again this can be understood if we take into 
account that the width of the conduction band decreases 
in this case and the electrons are transferred from the d 
to the / character orbitals (see [[]). In this way the sys- 
tem evolves from the itinerant to the localized situation 
(e/ = e a < Ep). According to our results (see Fig.[To|), 
this change should increase the value of the zero tem- 
perature magnetization and simultaneously decrease the 
Curie temperature T c {\AE\ is strongly reduced because 
the effective magnetic interaction in the localized limit, 
J rkky i is order V i [Q). This anomalous behavior has 
been experimentally observed by Shaheen et al |54| in 
La x Cei_ 2: Rh3B2. 

We can also connect our mechanism with the hydro- 
static pressure dependence of T c . To do this we calculated 
\AE\/N by the QMC method as function of increasing tb 
(Fig. [L4|a). Here we are assuming that the main effect of 
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the hydrostatic pressure is to increase % and to leave the 
other parameters unchanged. The order of magnitude 
of \AE\/N, which should be proportional to T c , and its 
qualitative behavior in Fig. 3a are in good agreement 
with the experimental results for CeRh 3 B 2 We see 
from Fig. 3a that for the itinerant FM case, |AU|/N is 
of the order of 100°K. This scale is much larger than the 
magnitude of the RKKY interaction |35| (~ 1°K) which 
is commonly used to explain the origin of the magnetic 
phase when the a electrons are localized. We also find 
that the FM state appears close to quarter filling and 
disappears for n close to 3/8. 
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FIG. 14. a) Energy difference per site between the FM and 
PM states as a function of tb- 



The other unusual property of these compounds is 
the shape of the magnetization curve versus tempera- 
ture which has a maximum below T c |66| . Again this is a 
property which can be easily explained (see the subsec- 
tion about the Ce compounds) within the context of the 
PAM. In addition, the order of magnitude of the Curie 
temperature of these compounds coincides with the en- 
ergy scale obtained from the PAM for the intermediate 
valence regime. 

The uranium monochalcogenides, like the Ce based 
compounds above described, exhibit a non-monotonic be- 
havior for the T c as a function of pressure |l9| . Fig. |l4| 
shows that this behavior can also be explained with the 
PAM. Notice however that the non-monotonic behavior 
shown in Fig. |lj has nothing to do with a competition 
between Kondo and RKKY interactions. 

Finally, the spin wave dispersions of these compounds 
also present some anomalies. For instance, neutron scat- 
tering experiments on a single-domain UTe crystal ( 
[^8|,^9|) show that for wave vectors q perpendicular to the 
ordered moment the excitations become more damped 
with increasing q. In US only a broad continuum of mag- 
netic response is observed |70[ ] . Damped and unpolarized 
spin waves are observed in USe J7l|-^3| . These properties 
indicate that the itinerant character of the / electrons is 
essential to have a good description of the magnetic ex- 
citations. 



C. Transition Metals 



B. Uranium Compounds 

The FM uranium monochalcogenides US, USe and 
UTe are semi-metals with large Curie temperatures of 
T c = 180A, 160A and 108A' and ordered moments of 1.5, 
2.0 and 2.2/ib respectively Most of the magnetic 

properties of these systems are still unexplained. The 
purpose of this subsection is to argue that the new mecha- 
nism for ferromagnetism introduced above is a good can- 
didate to explain some of the mysteries related to these 
compounds. 

Erdos and Robinson J6(| suggested that the uranium 
monochalcogenides are mixed valence systems. This sug- 
gestion was reinforced by measuring the Poisson ratio as 
a function of the chalcogen mass with low temperature 
ultrasonic studies on USe and UTe J67J . The coexistence 
of an intermediate valence regime and ferromagnetism is 
one of the unexplained properties of these compounds as 
it is recognized in Rcf . |24| . According to the traditional 
picture ]45| ] , the system should behave as a nonmagnetic 
collective Kondo state in the intermediate valence regime. 
In contrast to this picture, our results show that a par- 
tially saturated ferromagnetic state is stabilized in the 
mixed valence regime. This can explain the first striking 
property of the uranium monochalcogenides. 



Even though the transition metals are the most well 
studied itinerant ferromagnets, the ultimate reason for 
the stabilization of the FM phase is still unknown. Since 
the minimal correlated model (Hubbard Hamiltonian) 
proposed to describe these systems does not seem to have 
a FM solution, it is reasonable to ask whether an ex- 
tension of this minimal model, including more than one 
band, is necessary and enough to stabilize the FM solu- 
tion. The correlated 3d band of the transition metals is 
hybridized with weakly correlated and dispersive s and 
p bands. This situation is similar to the case already 
described for the / electron compounds. Therefore, it is 
natural to ask if there is a connection between the itin- 
erant ferromagnetism of the / and the d electron com- 
pounds. Notice that the order of magnitude of the T c 
is the same. Following the same motivation and using 
DMFT, Schwieger and Nolting Q concluded that the 
c?-band ferromagnetism can be stabilized when the hy- 
bridization between both bands is small. However, these 
authors also find a FM solution for the one band problem. 

In the case of the transition metals, the dispersion of 
the narrower band (3d) cannot be neglected. For this 
reason, we studied the stability of the FM solution as a 
function of t a . We can see from Fig. [l5] that the FM 
phase is even more stable for t a ~ —0.1 than for t a = 
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and becomes unstable for t a ~ 0.05. The reason for this 
asymmetric behavior is easy to understand in terms of 
the variation of S a : If t a is negative, then the effect of t a 
on the dispersion of the 4> band is opposite to that of the 
hybridization V (see Fig. [l6]). When t a ~ —0.1, we get, 
for the given e a and V, the minimum value for 5 a /A and 
therefore the most stable FM case. When we depart from 
this value of t a , S a /A increases, \AE\ decreases, and the 
FM state becomes less stable. 
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FIG. 15. Influence of the hopping t a on the FM state. 

This result indicates that the hybridization between 
bands can play a crucial role for the ferromagnetism of 
the iron group. In other words, the ferromagnetism of 
the transition metals can originate, at least in part, in 
the interplay between the correlations and the particular 
band structure and not solely in the intra-atomic Hund's 
exchange |iHj] . 
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FIG. 16. Band structure for different values of t a . 



The hybridization between bands and the particular band 
structure play a crucial role in this mechanism because 
they generate a multi-shell structure for the correlated 
orbitals. This structure, when combined with the a lo- 
cal Coulomb repulsion, favors a ferromagnetic state. The 
mechanism is analogous to the one which generates the 
atomic Hund interaction. In this sense, this is a gen- 
eralization to the solid of the atomic Hund's rule. The 
mechanism works in any finite dimension. 

The determination of a minimal model to explain the 
metallic ferromagnetism of highly correlated systems has 
been the object of intense effort during the last forty 
years. The results presented in this paper suggest that 
the PAM is a minimal Hamiltonian which can explain the 
itinerant ferromagnetism without including any explicit 
FM interaction. 

Another important aspect of this ferromagnetic solu- 
tion is its mixed valence character. According to the tra- 
ditional picture [ ft5| , the mixed valence regime should be 
a PM Kondo state. The appearance of a ferromagnetic 
instability in this region of doping rises some questions 
about the entire validity of a Kondo-like description in- 
spired by the one impurity problem. Even the PM phase 
obtained for U < U c is strongly influenced by the prox- 
imity to a FM instability (5§-||] . 

We have discussed the relevance of these results for 
some /-electron compounds which are itinerant fer- 
romagnets with high Curie temperatures (~ 100°K). 
In particular, the are several unusual characteristics 
of the Ce based compounds Ce(Rhi_a;Ru x )3B2 and 
La aj Cei_ 2; Rh3B2, and the uranium monochalcogenides 
US, USe and UTe, which can be explained, at least at 
a qualitative level, with the present mechanism. 

We have also considered the case relevant for the iron 
group where the dispersion of the lower band is not neg- 
ligible. The fact that the ferromagnetism is even more 
stable for finite values of t a when the hoppings of both 
bands have opposite signs indicates that our mechanism 
is relevant to explain the ferromagnetism of the transition 
metals, like Ni, where a correlated and narrow 3d band 
is hybridized with the 4s band. It suggests that the fer- 
romagnetism in the transition metals can originate, at 
least in part, in the interplay between the correlations 
and the particular band structure, and not solely in the 
intra-atomic Hund's exchange fig). 



VII. CONCLUSIONS 
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